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Abstract: We investigate the maximum size of CP violating effects in L>-mixing within 
the Standard Model (SM), using Heavy Quark Expansion (HQE) as theoretical working 
tool. For this purpose we determine the leading HQE contributions and also a s corrections 
as well as subleading l/m c corrections to the absorptive part of the mixing amplitude 
of neutral D mesons. It turns out that these contributions to Ti2 do not vanish in the 
exact SU(3)f limit. Moreover, while the leading HQE terms give a result for Ti2 orders of 
magnitude lower than the current experimental value, we do find a sizeable phase. In the 
literature it was suggested that higher order terms in the HQE might be much less affected 
by the severe GIM cancellations of the leading terms; it is even not excluded that these 
higher order terms can reproduce the experimental value of y. If such an enhancement is 
realized in nature, the phase discovered in the leading HQE terms can have a sizeable effect. 
Therefore, we think that statements like: "CP violating effects in D-mixing of the order of 
10~ 3 to 10 -2 are an unamhigous sign of new physics " — given our limited knowlegde of the 
SM prediction — are premature. Finally, we give an example of a new physics model that 
can enhance the leading HQE terms to Ti2 by one to two orders of magnitude. 
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1. Introduction 

Meson-antimeson oscillations have long since provided a rich area for theoretical studies. 
In an application of pure quantum mechanics, the concept of neutral kaon mixing [[jj led to 
the understanding of the observed CP violation in the decay K\, — > tttt |2|. Furthermore, 
hints to the mass of the c quark ||| were obtained from K mixing before the first evidence 
of the J/ty 0, Also Bd — B,i mixing H, Q provided information on the mass of the t 
quark prior to its discovery |8|, ^| . The latter was possible due to the sensitivity of meson- 
antimeson mixing to heavy virtual particles propagating in an internal loop of the transition. 
Even today, the absence of a Standard Model (SM) tree-level background turns precision 
measurements of meson-antimeson mixing into an excellent probe for new physics effects; 



recently, a possible indication of such effects in the phase <& s of the B s — B s system [10, 11 1 



has stirred vivid discussions |12, O, H| . 

Among the four mixing systems (K°, D°, Bj and B s ) the D° system is in a sense unique. 
The mixing mechanism relies on internal d-type quarks; due to the smaller mass of the b 
quark compared to the t quark, the kinematics of the dispersive part of the mixing amplitude 
are not completely dominated by the heavy third generation quark. Furthermore, due to 
the specific structure of the Cabibbo-Kobayashi-Maskawa couplings [jDj], the absorptive 
part Ti2 will feature an extremely efficient Glashow-Iliopoulos-Maiani (GIM) mechanism 



1 16]. We will discuss this in detail later on and show that it leads to a suppression of the 
leading contribution by several orders of magnitude. 

On the experimental side, the first evidence for charmed meson oscillation was reported 



by the Belle collaboration 0, |§, by BABAR |L9j and later by CDF [g0|. Currently, the 



relative decay width difference y and the relative mass difference x have been measured 
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to about 20% accuracy. The Heavy Flavor Averaging Group (HFAG) quotes the best-fit 
values pH | 

AT o AM , 9 r o . . 

V ■= = (7-3 ± 1.8) x 10- 3 , x := -— = 9.1± 2 | x lO" 3 . (1.1) 

The theory status for the D° system is, unfortunately, in a slightly worse shape. There 
are two main approaches to this issue: Heavy Quark Expansion (HQE) p2| , p3| , pi| is an 
expansion of the bilocal AC = 2 matrix elements as a series of local operators of increasing 
dimension, which are suppressed by powers of the heavy quark mass. While this technique 
provides an excellent tool to study B mixing |25||, the D° system predictions p6|, 27, 28, 29 1 



differ from experiment by up to a factor 10 3 . There are four main lines of argumentation 
why such a behavior is observed. First of all, a breakdown of the expansion in powers of the 
charm mass, which may not qualify it as "heavy" in the sense of the expansion, is possible. 
Comparing the ratios of the typical hadronic scale A and the heavy quark mass in the B 
and D° system (A/nib ~ 0.05 versus A/m c ~ 0.25), the expansion parameter has increased 
by almost a factor of five. Secondly, one can invoke a violation of quark-hadron duality 
due to non-perturbative long distance effects; reliable quark-level predictions for the decay 
width difference may consequently serve as a probe of quark-hadron duality in the charm 
system. Furthermore, such a deviation could arise, if the severe GIM cancellations present 
in the leading terms of the OPE are lifted for higher dimensional operators. Finally, new 
physics may enhance the SM result for D° — D° mixing. 

A second way to access the calculation of the D° — D° decay width difference Ar is based 



on exclusive techniques [30, 31 1. In principle one has therefore to determine all contributions 
to T D o and T^o with high precision, which is clearly beyond our current ability. As a first 
step to determine the size of Ar within the exclusive approach the authors of [30,31] take 
only the difference of the phase space of the corresponding final states into account. 

In this work we reexamine the HQE for the relative decay width difference y. We begin 
with a short introduction of the formalism for the D° — D° mixing via box diagrams. Sect. || 
deals with the leading-order (LO) HQE predictions for the mixing matrix element Ti2 in 
next-to-leading order (NLO) in a s . Furthermore, we will show that, contrary to expectation, 
the large cancellations due to GIM mechanism can generate a sizeable imaginary part in 
Ti2. After a brief discussion of possible effects of higher order terms in the HQE, a new 
physics model is presented, which can substantially enhance the leading HQE term. We 
finish with a conclusion. 



2. Mixing formalism 

The mixing of neutral mesons is described by box diagrams with the absorptive part Ti2 
and the dispersive part Mi 2 . The observable mass and decay rate differences are given by 
{<t> := arg[-M 12 /r 12 ]) 

(AM) 2 - J (Ar) 2 = 4|M 12 | 2 - |r 12 | 2 , 

AM Ar = 4|Mi 2 ||ri 2 | cos(0) . (2.1) 
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d,s d,s 

Figure 1: Contributions to T12 from operators of dimension 6 (D — 6). The leading order QCD 
diagram is shown in the left panel, an example for a s corrections is shown in the right panel. 



If |ri2/Mi2| <C 1, as in the case of the B s system (~ 5 • 10 3 ) or if (j) <C 1, one gets the 
famous approximate formulae 

AM = 2\M 12 \, Ar = 2|ri 2 |cos</>. (2.2) 

The experimental values for x and y suggest that in the D° system IT^/M^I ~ 0(1), the 
size of the mixing phase (j) will be discussed below. 



3. Leading HQE predictions 

The absorptive part of the box diagram with internal s and d quarks can be decomposed 
according to the CKM structure as 

Ti2 = — (^r ss + 2X s X d T sd + A^T rfrf ) , (3.1) 

with X x = V CX V* X . The application of the heavy quark expansion (HQE), which turned out 
to be very successful in the B system, to the charm system typically meets major doubts. 
Our strategy in this work is the following: instead of trying to clarify the convergence 
of the HQE in the charm system in advance, we simply start with the leading term and 
determine corrections to it. The size of these corrections will give us an estimate for the 
convergence of the HQE in the D° system. To this end we first investigate the contribution 
of dimension-6 (D = 6) operators to Ti2, see Fig. [l|. Next we include NLO-QCD corrections, 



which were calculated for the B s system [32, B3l pj, 13| and subleading terms in the HQE 



(dimension-7 operators), which were obtained in [35, 36]. To investigate the size of a s and 
l/m c corrections in more detail we decompose T ss s ddd the Wilson coefficients G and 
Gs of the AC = 2 operators Q and Qs (for more details see ]13|]) 

Tf = 6 ' 7 = [G-^IQI^ ) + GT (D°\Qs\D )] + ri, (3.2) 



with 



Q = u a ^(l - 7 5 )c a • 1^7^(1 - 75)0/3, 
Qs = u a (l + 75)c Q • ug(l + 75)173. (3.3) 
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The effect of the QCD corrections has already been discussed in [p^ j. In our numerics we 
carefully expand in a s : the leading order QCD contribution consists of leading order AC = 1 
Wilson coefficients inserted in the left diagram of Fig. (Eh, while our NLO result consists 
of NLO AC = 1 Wilson coefficients inserted in both diagrams of Fig. (|l|) and consistently 
throwing away all terms which are explicitly of 0(a 2 s ). Following ]38|| we have also summed 



terms like zlnz to all orders; therefore, we use in our numerics z = m s (m c ) /m c (m c 



0.0092. The matrix elements in Eq. (3.3) are parameterized as 
(D°\Q\D°) = lf D oM 2 DO B^), 



(3.4) 



(^IQsl^ ) 



f D oM 



M D o 



mjji) + m u (n) 



We take f D = 212(14) MeV from |9) and we derive B(m c ) = 0.9 ± 0.1 and B s (m c ) ~ 1.3 
from p0| , [LT[| . Finally we use the MS scheme for the charm mass, m c (fh c ) = 1.27 GeV. For 
clarity we show in the following only the results for the central values of our parameters, 
the error estimates will be presented at the end of the next section. We obtain 
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For the error estimate we vary fi\ between 1 GeV and. 2?7i c . Combining G and. G$ to r a 
we get (in units of ps _1 ) 





LO 


ANLO-QCD 


Al/m c 






3.52 


-0.94 (-27%) 


-0.76 (-22%) 


1.82 




3.54 
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1.87 



Using instead the operator basis suggested in |L3[ with 

Qs =u a (l + 75)0/3 • up{l + 75)c Q , 



(3.5) 



(D°\Qs\D°) 

and Bs(m c ) ~ 1.2 leads to 



M D o 



m c (n) + m u (n) 
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All in all we get large QCD (up to 50%) and large l/m c corrections (up to 30%) to the 
leading D = 6 term, which considerably lower the LO values. In the (Q, Qs)-basis numercial 
cancellations can occur which mimic very small QCD corrections for T xy . Despite large 
corrections, the HQE seems not to be completely off. From our above investigations we see 
no hints for a breakdown of OPE. The same argument can be obtained from the comparison 
of B and D meson lifetimes. In the HQE one obtains 

tb = r o,po + 5Tp° ^ r o,D" L | ST D o \ A 5T B \ 

where the leading term Tq oc mjj c Vq KM corresponds to the free quark decay and all higher 
terms in the HQE are comprised in 5T. For the ratio r jjo/Tq^b one gets a value close to 
one. Higher order HQE corrections in the B system are known to be smaller than 10 % 



1 25]. Using the experimental values for the lifetimes we get 



_ 1.4...4(Exp.)«l- 1 + -^ . (3.7) 

From this rough estimate one expects higher order HQE corrections in the D system of up 
to 300 %. So clearly no precision determination will be possible within the HQE, but the 
estimates should still be within the right order of magnitude. 



4. Cancellations 

As is well known huge GIM cancellations [ |l6| arise in the leading HQE terms for D° 
mixing. To make these effects more obvious, we use the unitarity of the CKM matrix 
(Xd + X s + Xb = 0) to rewrite the expression for the absorptive part in Eq. (|3.1| ) as 

T 12 = -X 2 S (T ss - 2T sd + T dd ) + 2X s X b (T sd - T dd ) - X 2 b T dd . (4.1) 

Note that the CKM structures differ enormously in their numerical values: A^ jS oc A and 
Xb oc A 5 in terms of the Wolfenstein parameter A ~ 0.2255. In the limit of exact SU(3)f 
symmetry, T ss = T sd = T dd holds and therefore, contrary to many statements 1 in the 
literature, Ti2 = — A^r^ is not zero although strongly CKM suppressed. Next we expand 
the arising terms in z. Using Eq. ( [3.2; ) we get in LO 

r£= 6 = 3.55477 - 3.22581z - 14.877z 2 + ... , 

rf d = 6 = 3.55477 - 1.612911 - 7.43849z 2 + ... . (4.2) 

The first term in the above equations obviously corresponds to r^ =6 . For the combinations 
in Eq. ( [4.1[ ) we get 

(T ss - 2T sd + T dd ) D=6 = -36.91z 3 « A 7 ' , 

(T sd - T dd ) D = e = -1.6131 « A 28 . (4.3) 

1 These statements are obtained assuming V u b = 0, which we have shown to be not justified for the HQE 
approach. 
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To make the comparison with the arising CKM structures more obvious, we have expressed 
the size of these combinations also in terms of powers of the Wolfenstein parameter A. As 



is well known, we find in the first term of Eq. (4.1) an extremely effective GIM cancellation, 
only terms of order z 3 survive. In NLO we get 



r 



1.8696 



T D ~ Q ' 7 = 1.8696 



5.5231z 
2.7616z 



13.8143z 2 + 
7.4906z 2 + .... 



(4.4) 



The arising combinations in Eq. ( f4.1| ) read now 

^= 6 ' 7 = 1.17z 2 - 59.5z 3 « A 6 - 2 - A 6 - 7 , 
z '= 6 ' 7 = -2.76z«A 2 - 5 . 



(r ss — 2T sd + T dd ) 



(r sd -T dd r=°'< = -2.76z«A 2 - 5 . (4.5) 

The fact that now the first term of Eq. (f4.l| ) is of order z 2 compared to z 3 in the case of 



the LO-QCD value was discussed in detail in |33| and later on confirmed in [37|. These 
numbers are now combined with CKM structures, whose exact values read 



-C12C23C13S12 - cl 2 ci3S23Si3e tSvi = O (A 1 + zA 5 ) , 



+C12C23C13S12 - Sl2Cl3S23Sl3e i<513 

ci 3 s 23 s 13 e i5 ^ = O (A 5 + iA 5 ) , 



O (A 1 + i\ 7 ) , 



(4.6) 



with Cij = cos(9ij) and = sin(0y). Looking at Eq. ( [4.6D , it is of course tempting to 
throw away the small imaginary parts of A^ and X s , but we will show below that this is not 
justified. Doing so and keeping only the leading term in the CKM structure (C12C23C13S12), 
which is equivalent to approximate A& = 0, one gets a real Ti2 which vanishes in the exact 
SU(3)f limit. Keeping the exact expressions, we see that the first term in Eq. ( |4.l| ) is 
leading in CKM (O [A 2 + iA 8 ] ) and has a negligible imaginary part, but it is suppressed by 
A 6 ' 2 . The second term in Eq. (hi) is subleading in CKM (O [A 6 + iA 6 ] and it can 



have a sizeable phase. This term is less suppressed by SU(3)f breaking (~ 2.1 z ~ A 2 ' 5 ). 
The third term in Eq. ( |4.1[) is not suppressed at all by SU(3)f breaking, but it is strongly 
CKM suppressed (C [A 10 + iA 10 ]). For clarity we compare the different contributions of 



Eq. (4.1) in the following table 





1 st term 


2 nd term 


3 rd term 


LO 


37A 2 z 3 


3.2A s AbZ 


3.6A 2 




«A 90 


pa A 8 - 4 


pa A 9 - 1 


NLO 


1.17A 2 z 2 


5.5A s AfeZ 


1.87A 2 




«A 8 - 2 


« A 8 - 


pa A 9 - 6 



From this simple power counting, we see that a priori no contribution to Eq. (^~J) can 
be neglected. Taking into account the hierarchy of the CKM matrix elements 3 , we find 



2 With 'NLO' we denote the sum of leading order value, 0(a s ) QCD corrections and subleading — 
contributions. 

3 Actually \V u b\ is numerically of order A 4 and therefore A(, oc A 6 (see but in the literature V u b = 

AA 3 (p — irf) with small values of p and r\ is commonly used. 
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that the first two terms of Eq. (4.1) are of similar size, while the third term is suppressed. 
Moreover, the second term can give rise to a large phase in Ti2, while the first term has 
only a negligible phase. To make our arguments more solid we perform the full numerics 



using the CKM values from |43| and obtain for the three contributions of Eq. ( |4.1[ ) 

10 7 rf 2 =6,r = - 14.6409 + 0.0009i (1 st term) 

- 6.68 - 15. 8i (2 nd term) 

+ 0.27 - 0.28i (3 rd term) 

= -21.1- 16.0i = (11...39)e- i(a5 - 2 - 6) . (4.7) 

Here we show for the first time the errors; they are estimated by varying Hi between 1 GeV 
and 2m c and by taking into account the results for both choices of the operator basis. The 



first term in Eq. (4.1) turns out to be very sensitive with respect to the the exact values 
of the bag parameters and its real part is approximately of the same size as the second 
term, which features a large imaginary part. Furthermore, even the third term can give a 
non-negligible contribution, in particular to the imaginary part. 

To summarize, we have demonstrated that the typical approximation ~ 0, which 
is equivalent to neglecting the imaginary parts of Ad and A s is wrong for the case of the 
leading (D=6,7) HQE prediction for y and yields the wrong conclusion that p^ _6 ' 7NLO 
cannot have a sizeable phase. We get for the first terms in the OPE a value for y of 

y D=6,7NLO < |p i2 | . ^ = 4 7 . 10 -7 L g . 1Q -6 _ ( 4 g) 

The range of values was again estimated on the basis of the renormalization scheme depen- 
dence and the choice of the operator basis. These values are still a factor of 0.5. ..1.6 • 10 4 
smaller than the experimental number. This is in contrast to our previous expectations that 
the HQE should give at least the right order of magnitude. Moreover, we do not confirm 



the observation made in [37| that the NLO result for is almost an order of magnitude 
larger than the LO result. 



5. Higher HQE predictions 

higher order terms in the HQE of D° mixing were discussed. If the GIM 
cancellation is not as effective as in the leading HQE term, operators of dimension 9 and 
dimension 12, see Fig. (|2|), might be numerically dominant. In order to obtain an imaginary 
part of the loop integral, the operators of dimension 9 have to be dressed with at least one 
gluon and the operators of dimension 12 with at least two gluons. If we normalize the 
leading term (left figure of Fig. (|l])) to 1, we expect the D = 9 diagram of Fig. (|2[) to be of 
the order 0(a s (47r)(qq) /m^) ~ 0.03 and the D = 12 diagram of Fig. (§) to be of the order 
0(a1(4TT) 2 (qq) 2 /m®) ~ 10~ 3 . As explained above, the formally leading term of D = 6 is 
strongly GIM suppressed to a value of about 2 • 10 -5 and the big question is now how severe 
are the GIM cancellations in the D = 9, 12 contribution. For the contributions to y we get 



In pq, 27, 28 
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{(id), (ss) 



{dd), (ss) 




d,s (dd),(ss) 

Figure 2: Contributions to Ti2 from operators of dimension 9 (D = 9, left panel) and dimension 
12 (D = 12, right panel). To obtain an imaginary part the D — 9 diagrams have to be dressed with 
at least one gluon and the D — 12 diagrams with at least 2 gluons. 
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If there would be no GIM cancellations in the higher OPE terms, then the D = 9 or D = 12 
contributions could be orders of magnitudes larger than the D = 6 term, but in order to 
explain the experimental number still an additional numerical enhancement factor of about 
15 has to be present. For more substantiated statements T^ -9,12 has to be determined 
explicitly, which is beyond the scope of this work [[|4) . This calculation is also necessary in 
order to clarify to what extent the large phase in T\2 from the first OPE term will survive. 
In order to determine the possible SM ranges of the physical phase 0, in addition one has 
to determine Myi- 



6. New physics 

Finally we would like to address the question, whether new physics (NP) can enhance T^. 
In the B s system it is argued J45|| that T\2 is due to real intermediate states, so one cannot 
have sizeable NP contributions. Moreover, the mixing phase in the B s is close to zero, 
so the cosine in Eq. ( |2.2j ) is close to one and therefore NP can at most modify 0, which 
results in lowering the value of AT compared to the SM prediction. In principle there is a 
loophole in the above argument. To T\2 also AB = 1 penguin operators contribute, whose 
Wilson coefficients might be modified by NP effects. But these effects would also change 
all tree-level B decays. Since this is not observed at a significant scale, it is safe to say 



that within the hadronic uncertainties Ti2 = Tf^ 1 and therefore the argument of |45| holds. 
Since in the D° system the QCD uncertainties are much larger, also the possible effects 
might be larger but not dramatic. The peculiarity of the D° system - the leading term in 
the HQE is strongly suppressed due to GIM cancellation - gives us however a possibility 
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5 10 15 20 25 30 35 



Figure 3: The enhancement factor rf^ 14 /]?^ 13 using the possible values for Vckm4 found in | |42| . 
The color encoded scale denotes the logarithm of the number of allowed parameter points of Vckm4- 



to enhance Ti2 by a large factor, if we manage to soften the GIM cancellation. This might 
be accommodated either by weakening the SU(3)p suppression in the first two terms of 



Eq. (4.1), see, e.g., Petrov et al. [|29|], or by enhancing the CKM factors of the last two 
terms in Eq. (f4.1| ). The latter can be realized in a model with an additional fourth fermion 
family (SM4). The usual CKM matrix is replaced by a four dimensional one (VckM4) and 
the unitary condition now reads + A s + + Ay = 0. Eq. ([4. 1|) is replaced by 



12 — — A^ (T ss — 2T sd + T a 
+2A s (Af, + Ay) (T S( j - 



3.1) 



dd ) 



(Ah + Afe/) 2 r a 



In p2|| an exploratory study of the allowed parameter space of VckM4 was performed and 
as expected only very small modifications of A^ and A s are experimentally allowed. In 
almost all physical applications these modifications are numerically much smaller than the 
corresponding hadronic uncertainties and therefore invisible. However, in the D° mixing 
system it might happen, that all dominant contributions cancel and only these modifications 
survive. In the SM the first two terms of Eq. fl4.l|) are numerical equal. In the SM4 the 
numerical hierarchy depends on the possible size of Ay, see Eq. fl6.ll) . In particular, it was 
found in |f42|| that currently a value of Ay of the order A 3 is not excluded. This means that 



the second term of Eq.(6.1) could be greatly enhanced by the existence of a fourth family 
and also the third term would now become relevant. Using experimentally allowed data 
points for VckM4 from |f42[ we have determined the possible values of Ti2 in the SM4 from 
Eq. ( |6.1P ; enhancement factors of a few tens are possible, see Fig.(|3|). It should be noted 
that this enhancement is very sensitive to the exact values of the bag parameters; e.g. a 
reduction of B$ by 30% triples the range of possible enhancement factors. 



7. Discussion and Conclusions 

In this work we investigated the leading HQE contribution to the absorptive part of D 
mixing and the leading corrections to it. We found that the size of these corrections is 



- 9 - 



large, but not dramatic (~ 50% QCD, ~ 30% l/m c ). So we see no signal for a breakdown 
of the OPE and it seems that the HQE might be appropriate to estimate the order of 
magnitude of T\2- For a further investigation of the question of the convergence of the 
HQE for the case of the D system a systematic study of D meson lifetimes within that 
framework might be very helpful. 

Above we have explained in detail that r{^ — 6 ' 7 gives, due to huge GIM cancellations, a 
value of y which is about a factor of about 10000 smaller than the current experimental 
expectation, but it can have a large phase and it also does not vanish in the exact SU(3)f 
limit. The possibility of a sizeable phase in T\2 is a new result. The important question is 
now, how big can the physical phase <p be? 

Due to the peculiarity of the D° system — the extreme GIM cancellations — it might also be 
possible that the HQE result is dominated by D = 9 and D = 12 contributions, if there 
the GIM cancellations are less pronounced, see e.g. To quantify that possibility these 

higher dimensional corrections have to be determined explicitly, i.e beyond the estimates 
presented in [27], This calculation will also show what values are possible for the 
phase of T^. To determine finally the physical phase one has to determine in addition 
M 12 . 

Currently, estimates about the possible size of the phase in D mixing are typically 
based on the assumption that the A^-term can be neglected. The remaining contribution 
proportional to A^(r ss — 2T S d + T dd leaves almost no space for a CP-violating phase. In this 
paper we have shown that this widely used assumption is wrong for the leading (D=6,7) 
HQE term: the A^-term is as sizeable as the pure \ s term and it can have a huge phase. 
Although we can not proof at the current stage that a sizeable phase can survive after- 
all corrections are included in the HQE calculation, we think one should meet claims that 
already small values for a D mixing phase are an unambiguous sign for new physics with 
some caution. This question has to be studied in more detail both within the inclusive 
approach — as explained above — and within the exclusive approach. 



To become more concrete, let us speculate: if the first term in Eq. (4.1) is enhanced 
by a factor r/ 3 due to higher dimensional terms in the HQE, where the GIM-suppression 



is much weaker compared to the leading HQE terms, as advocated e.g. in [28|, the same 
effect — but less pronounced — will also be active in the second term of Eq. (4T). For the 
first term, the authors of [ f28| see a possibility of enhancing the effect of order z 2 up to \fz 
(by 3 powers of \fz). For the second term we expect an enhancement from z to ^fz (by 1 
power of ^fz) . So the enhancement of the second term of Eq. ( |4.1[) is a factor of the order 
of r]. Now we can make some numerical predictions for y and Im T^/ Re Ti2 depending on 
the numerical enhancement factor r\ - we use only values for r] that are within the estimates 



V 


lmr 12 
Rer^ 


^Theory 


Deviation from exp. central value of y 


18.24 


0.32% 


0.73% 


0a 


14.54 


0.50% 


0.37% 


2a 


11.64 


0.78% 


0.19% 


3a 


4.33 


5.62% 


0.01% 


4a 



(7.1) 
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We see that the new contribution, we have found in this paper could lead to a relative size 
of the imaginary part of Ti2 of the order of up to 1 %, while the absolute value of T12 is close 
to the experimental value of y. A size of the physical phase of the order of 1% within the 
standard model is typically excluded in the literature, where the correction we have worked 
out in this paper is not taken into account. But we would like to warn the reader: for the 
numbers in Eq. ([7T]) we have purely speculated about the size of the enhancement factor 
7] - but they are within the estimates of ]28| |. Its concrete value has to be determined by a 
calculation of the higher order terms in the HQE. This task, however, is beyond the scope 
of the current paper; we plan to investigate this question in [44|. 

Finally we have shown that new physics, in particular a small violation of the unitarity 
of the 3x3 CKM matrix, can enhance the leading HQE prediction for Ti2 by a double-digit 
factor. 
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